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Abstract: We derive the evolution of a hadronic hght cone wave function with energy at weak 
couphng. Our derivation is vahd both in the high and the low partonic density limit, and thus 
encompasses both the JIMWLK and the KLWMIJ evolution. The hadronic wave function is shown 
to evolve by the action of the Bogoliubov-type operator, which diagonalizes on the soft gluon sector 
the light-cone hamiltonian in the presence of an arbitrary valence charge density. We find explicitly 
the action of this operator on the soft as well as the valence degrees of freedom of the theory. 
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1. Introduction. 

The problem of calculating hadronic scattering amplitudes at high energy is an old one. It goes 
back to classical works of Gribov on Reggeon Field Theory [0 in the pre-QCD days. Within the 
framework of QCD this question has been addressed from different points of view p|, ^, |^ . 

In the last ten years or so the subject has seen new developments. Some of these have been 
triggered by Mueller's reformulation of the BFKL equation Q in terms of the dipole model 
with additional input provided by the functional approach of [^. The result was the derivation 
of the functional evolution equation for the hadronic amplitude - the so-called Balitsky hierarchy 



1^ or JIMWLK equation |jl0|, |Tl|; and its simplified mean field version due to Kovchegov ||12|| . 
This evolution takes into account coherent emission effects in the dense hadronic wave function, 
or in other words partonic saturation effects. These effects lead to unitarization of the scattering 
amplitude. Although the language of this approach is different from the original Reggeon Field 
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Theory, a direct relation between the JIMWLK evolution and the QCD Reggeon Field Theory has 
been investigated recently []T3|, 



In the last couple of years spurred on by observations of |T5[, the realization has emerged 
that the existing evolution equations, which are tailored to describe the situation when a small 
perturbative projectile scatters off a large dense target, do not include the so called Pomeron loop 
effects. The effort to account for the Pomeron loops using the probabilistic view of the evolution |]TB| 
has lead to interesting analogies between QCD and statistical systems Alternative approaches 
based on effective Lagrangian can be found in [0, |0|, |li 



Another avenue that has been explored in this context is the direct approach to the evolution 
of the hadronic wave function [^, ^ ^ This approach yielded the evolution equation valid 
in the limit opposite to that of JIMWLK, namely when the hadronic wave function is dilute. 
This so called KLWMIJ equation pi| is related to the JIMWLK equation by the dense-dilute 



duality transformation The basic strategy of this approach is to calculate the light cone 

hadronic wave function of soft gluonic modes, given the color charge density j°'{x) due to the 
'valence' modes - the modes with large longitudinal momentum. When the hadron is boosted, 
the longitudinal momentum of the soft modes is increased and they contribute to the scattering 



matrix and other physical observables [27|. The evolution of any physical observable is therefore in 
principle completely determined once we know the soft part of the wave function. 

So far the hadronic wave function has been calculated only in the KLWMIJ limit, namely 
when the valence charge density is small; ~ g. The JIMWLK evolution on the other hand 

is valid when j°'{x) ~ ^/g, but no wave function evolution is available in the JIMWLK regime. 



The derivation of is given directly for the scattering matrix. The original derivation of ||10 
is not far in spirit from the wave function form of the evolution, however it involves additional 
approximations which do not allow to read off the evolution of the wave function directly from the 
JIMWLK equation. 

The main motivation to know the explicit form of the wave function evolution comes from the 
possibility to use it to derive the generalization of the JIMWLK/KLWMIJ evolution that includes 
the Pomeron loops. The knowledge of the wave function is also crucial to be able to address a wide 
range of semi-inclusive observables p5|. 

In the present paper we derive the soft gluon wave function valid at any physically interesting 
value of the valence color charge density. The expression we derive is valid both in the JIWMLK 
and the KLWMIJ limits as well as at any value of the valence charge density which interpolates 
between the two: g < j < 1/g . We do this by diagonalizing the leading part of the light cone 
Hamiltonian on the soft gluon sector. The transformation that diagonalizes the Hamiltonian turns 
out to be of the Bogoliubov type with parameters depending on the valence color charge density 
operator. We find explicitly the action of this transformation on the basic quantum degrees of 
freedom: the soft components of the vector potential A^{x, x~) and the valence color charge density 
fix). 

We show that the expression for the wave function indeed reproduces the JIMWLK and the 
KLWMIJ evolution equations. To reproduce the KLWMIJ equation one simply neglects the coherent 
emission effects in the wave function. Thus the evolution of the wave function in this limit is strictly 



perturbative 1 24]. The nontrivial physics in this limit is entirely due to the multiple scattering 
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corrections in the scattering amplitude. On the other hand to derive the JIMWLK hmit we keep 
all the coherent emission effects in the wave function. However as we show explicitly below, in this 
limit we neglect certain multiple scattering corrections in the scattering amplitude. Physically this 
is justified in the situation where our hadron scatters on a perturbatively small target, which is 
when the JIMWLK evolution is valid. 

To derive the evolution equation which includes Pomeron loops exactly we have to keep both 
types of effects in the evolution of the scattering amplitude. Within the present framework this 
looks like a tractable problem. It is however beyond the scope of the present paper and is left for 
future work. 

The paper is structured as follows. In Sec. II we recall the general framework of the high 
energy evolution. Sees. Ill and IV are the main part of this paper. Sec. Ill is devoted to the 
derivation of the "vacuum" wave function of the soft gluon Hilbert space in the presence of the 
valence color charge density. In Sec. IV we show that this diagonalization is achieved by the action 
on the free vacuum of a Bogoliubov type operator and derive explicitly the action of this operator 
on the soft and valence degrees of freedom. In Sec.V we show how both the JIMWLK and the 
KLWMIJ evolution equations follow from the wave function we have found in Sec. Ill in different 
limits. Finally a discussion is presented in Sec. VI. 



2. High energy evolution 



The logic of our approach is the same as described in ||2^, Suppose that at some initial rapidity 
Yq we know the wave function of a hadron. In the gluon Fock space it has a generic form (we work 
in the A" = gauge) 

\P)y, = na^%x,k^)]\0) . (2.1) 

There is some minimal longitudinal momentum = A below which there are no gluons in this wave 
function. More precisely, the number of soft gluons with fc"*" < A is not zero but is perturbatively 
small so that their contribution to the scattering amplitude at Yq is a small perturbative correction 
and can be neglected. 

We are interested in describing the scattering of this hadron on some target. The target is 
described by some distribution of color fields ax = with a probability density distribution 
Wr[aT]. The second- quantized 5*- matrix operator in the eikonal approximation (in which we are 
working throughout this paper) is given by 



where 



S 



fix) = g 



exp 



i / (fxf{x)a^{x) 



a^\x, k+)T;^^a%x, k^ 



dk^ 



(2.2) 



(2.3) 



'fc+>A 27r 

is the color charge density operator at the transverse position x ( with = if"-'^'^ - the generator 
of the color group in the adjoint representation). After scattering on a particular configuration of 
the target field the hadronic wave function becomes 



S\P)y, = nS''\x)a^\x,k+)]\Q) 



(2.4) 
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where is a unitary matrix - the single gluon scattering matrix. Since the scattering amphtude 

is sensitive only to the color charge density in the hadronic wave function and not to any other 
characteristic of the hadron, we can think of this wave function as being specified by some distri- 
bution of j°'{x). The color charge density correlators are determined in terms of the 'probability 
density functional' W[j] via 



The forward scattering amplitude is then given by 



where 



S^Jar] = m*[a{x,k+)]Sna^{x,k+)]\0)= J Dj WyM exp 



i / d xf{x) 



(2.5) 



(2.6) 



(2.7) 



The total rapidity of the process is Y while the target is assumed to be evolved to rapidity Y — Yq. 
Here, characterizes the distribution of color fields in the target, while W [j] characterizes 
the distribution of color charges in the projectile. Due to Lorentz invariance S is l^o independent. 
The evolution of the S'-matrix ( p.6|) with energy in the high energy limit has the generic form 



d 



dY 



S = j Da''TW^_y^[aT{x)] H^^^^ 



6 



6 ax 



(2.8) 



where H^^^ is the Hermitian kernel of high energy evolution, which can be viewed as acting either 
to the right or to the left: 



_d_ 

"dY 



H 



RFT 



6 



6ai 



— 

dY 



H 



RFT 



6 



6ai 



W^lar] . (2.9) 



The color charge density operators are the generators of the SU{Nc) algebra and as such do 
not commute 

[fix], fiy)] = tr'^fix)6'ix-y). 

As explained in detail in [^, to properly take into account the non commuting nature of the 
charge density operators j{x) and to still be able to represent wave function averages in terms of 
the functional integral over 'classical' fields j", one has to assign to j an additional 'longitudinal' 
coordinate. Thus in effect j'^{x) — > j"'{x, x~), where the value of x~ simply keeps track of the order 
of the operators j in the correlation function eq.( |2.5| ). An analogous 'longitudinal coordinate' should 
be assigned to the target field a^- Since in this paper we work in the Hamiltonian formalism and 
explicitly keep track of the commutation relations of the quantum operators j°'{x), we will not need 
to dwell on this additional longitudinal coordinate. 

The preceding discussion is given in the situation when the increase of rapidity is assigned to 
the target. One can equally well boost the projectile. The evolution of the projectile probability 
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density functional W[j] is related to that of since the two are related by the functional Fourier 

transform eq.( p2.7D 



S_ 



-J 



W[j] . 



(2.10) 



As the hadron is boosted by rapidity AY, the longitudinal momenta of the gluons in its wave 
function are scaled by the boost parameter fc"*" — >■ e^^k~^. Thus some gluons in the wave function 
emerge after boost with the longitudinal momenta above the cutoff A and have to be taken into 
account in the calculation of the scattering amplitude. The number of thus 'produced' additional 
gluons in the wave function is proportional to the total longitudinal phase space / = AY. 

To find the evolution of the scattering amplitude we need two ingredients. First we have to 
solve for the initial hadronic wave function with greater accuracy on the soft gluon Hilbert space 
than is necessary to calculate the scattering amplitude at the initial rapidity Yq. Second we need to 
take into account the contribution of these soft gluons into the scattering amplitude at the rapidity 
Y = Yq + AY, which amounts to the transformation 

r-A 



fix) 



+ J 



soft 



J 



soft 



X 
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(2.11) 



in eq.( p.2|) . This transformation is conveniently represented in terms of the charge density shift 
operator (which also has the meaning of the 'dual' to the Wilson line operator] 



Ra 



exp 



d^zfsoftiz) 



6 



6fiz) 



fix) 



Rafix) 



(2.12) 



The crucial part of this program is the knowledge of the wave function on the soft gluon part 
of the Hilbert space, A;"*" < A with some minimal accuracy. The calculation of this wave function is 
the subject of the next section. 

The QCD light cone Hamiltonian H responsible for the dynamics of the soft modes is diago- 
nalized by the action of a unitary operator Qay, where AY corresponds to the phase space volume 
occupied by the soft modes. Equivalently, the vacuum wave function of the soft modes in the pres- 
ence of the valence color charges is Q\P). The kernel of the high energy evolution, H^^^ is related 
to VL as |l^ : 

(Oa|l^iy(j,a) (ka - l) ^^Ay (j, «) |0a) 

HRFt^_ 1 (2.13) 

Ay-.o AY ^ ' 

We will find below that VL is an operator of the Bogoliubov type for any physically interesting j: 

Vl = CB 

with C denoting a coherent operator, which is the exponential of an operator linear in the soft fields 
A, whereas B is an exponential of an operator quadratic in A. In the dilute limit j ^ g we have 
B = 1 and the coherent operator C leads to the KLWMIJ evolutionp^. For dense systems j ~ 1/g, 



the Bogoliubov operator B also contributes to the leading order evolution kernel H^^^ . We derive 



the action of on both the valence and soft degrees of freedom, which enter equation ( p.l3|) . The 
JIMWLK Hamiltonian |Ty] is obtained from the general expression ( p.l3| ) in the limit of weak target 
fields ax expanding Ra to second order in S/Sj. 
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3. Diagonalizing The Soft Gluon Hamiltonian 



We will proceed in the following steps. In section |3.1| we formulate the light-cone Hamiltonian for 
soft gluon modes fc+ < A, coupled to the color charge density of the hard modes > A. We 
observe that the zero modes of the vector potential are not independent degrees of freedom, but are 
constraint by the residual gauge fixing and the requirement of finiteness of energy. In section |3.2| , 
we solve the resulting constraints. In section p.3| , we diagonalize the resulting Hamiltonian, by first 
finding the complete set of solutions to the classical equations of motion, and then expanding the 
field operators in this basis. To ensure the canonical commutation relations for the creation and 
annihilation operators associated with these basis functions, a proper normalization of the classical 
solutions is needed. This normalization is found in section 



3.1 The Hamiltonian and the canonical structure 

The starting point of our approach is the light cone hamiltonian of QCD 



H = [ ^d'x (l-U-{k^,x)U-{-k+,x) + ]G^^{k+,x)G^^{-k+,x)) , (3.1) 
Jk+>o 271 \2 4 J 

where the electric and magnetic pieces have the form 
n-ix-,x) = -^{D'd+A,yix-,x), 

Gr(x-, x) = d^At{x-,x) - 9M;:(x-, x) - gr"^ 4(x-, x) Al{x-,x) , (3.2) 

and the covariant derivative is defined as 

Df 'i>^ = {di6^^ - gr^Af) <^>^ (3.3) 

Our aim is to diagonalize this Hamiltonian on the Hilbert space of soft gluon modes - those 
with longitudinal momenta smaller than some scale A. We assume that the valence part of the wave 
function (the component of the full wave function which does not contain soft modes) is known and 
is completely specified by the correlation function of the color charge density 

fix) ^ zgr"^ / -- a\\k\ X) atik^, x) . (3.4) 

The soft modes are the interesting dynamical degrees of freedom of our problem, and they interact 
with the valence ones via eikonal coupling in the Hamiltonian. The Hamiltonian for the soft modes 
is then given by eq.( |3.1| ) with the substitution 



+9 r'" I ^^lik"- - p\ x) (-ZP+) 4^(p+, x) . (3.5) 

-i{k+ + le) J\p+\<A 2tt 
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The soft fields A are defined only below the longitudinal momentum cutoff A, but we will not 
explicitly indicate it in the following. 

The canonical structure of the theory is determined by the commutation relations of the fields. 
As we will see, the zero momentum mode of the field A is non dynamical and is determined by 
the residual gauge fixing (still not specified so far on top of the usual light cone gauge condition 
A'^ = 0) and the constraint of finiteness of energy. We denote by A the part of the field that does 
not contain the mode with vanishing longitudinal momentum - the zero mode. 

The canonical commutators of the field A are 



(3.6) 



with 



One defines the light cone canonical creation and annihilation operators as usual through 



(3.7) 



A; (x . x) 



dk+ 1 



/o 27r y/2k+ 

'at{k+,x),ay{p^,y)] = {27r) 6^' 6,, 6{k+ - p^) 6'^'\x - y) 
This translates into (fc+ ^ 0): 



(3.8) 



At{k-',x),A%p+,y) 



TT 



1 1 

+ 



2 \k+ + ie k+ - it 



5{k+ + p+) 5""^ 5ij 5^^\x - y) 



(3.9) 



The Hamiltonian eq.( p.l|) commutes with the generator of the x - independent gauge trans- 
formation, which on physical states should vanish: 



j dx-iDid+A,r-f{x) = 0. 



(3.10) 



Following the standard procedure we should fix this residual gauge freedom by imposing a gauge 
fixing condition. We will be working in the gauge (same as in [|T0|] ) 



d.A'^ix- 



-oo 



0. 



(3.11) 



/^From previous analysis of the behavior of the field in this gauge |10|, we know that the vector 
potential vanishes at x" — oo but approaches a non vanishing asymptotic value at x~ — >■ oo, 
which we denote by A^{x^ oo,x) = 7f(x). Separating the nonzero momentum modes, we thus 
write 

A^{x-,x) = ljtix) + A^ix-,x). (3.12) 



-7- 



Even though A has no zero momentum mode, its asymptotics is not vanishing but is rather given 
by ±|7f • It is thus convenient to define a field c which has regular behavior at infinity by 

A'l{x-,x) = e(x-)7f(x) +c^(x-,x), 

c^(a;- ^ ±oo,x) ^ 0, (3.13) 

and 

d+A'^ = = 5ix-)-ft + a+c^ . (3.14) 

Our aim is to find the ground state of the Hamiltonian eq.( p.l| ) given the charge density j"" (more 
precisely we consider the matrix elements of the operators j"(x) on the Hilbert space of the valence 
modes as known). 

Our first observation is that since the Hamiltonian is the integral of the positive definite Hamil- 
tonian density over x~ , the necessary condition for finiteness of energy is vanishing of the density 
at x~ — > ±oo. The finiteness of the magnetic part of the Hamiltonian requires 

G^/x-^oo) = 0, (3.15) 



while the finiteness of the electric part is ensured by eq. (|3.10|) . 

We will use the gauge fixing condition and the finite energy conditions as operatorial constraints 
that determine 7 in terms of A (or equivalently q). This is equivalent to Dirac bracket quantization 
of the fields A which leave the canonical commutators of A unchanged. The commutators of 7" 
with Af and between themselves are then determined by solving the constraints. 

Expressing the magnetic constraint equation ( ^.15] ) in terms of 7", we obtain 

^]{x) - dj 7f (x) - gr"' i^{x) 7;(x) = . (3.16) 
To express the electric constraint eq.( p.lO| ) we use the fact that given the boundary conditions on 

dx-d+A^ = jf. (3.17) 
We then find^ 

- \gr'Hit - ar'^ j dx- A\{x-)d+ A'lix-) = fix) , (3.18) 

or, equivalently, 

d.^tix) - \gr'H{^m^)+gr''H{x),4{xM~9r''' j dx-c\{x-)d+ct{x-)=3\x). (3.19) 

In this equation c(0) should be understood as 

ct{x~ = 0) = A{x- = 0) = ^K(x- = 0+) + c1{x- = 0-)] , (3.20) 

where 0"*" = + e; 0~ = — e; e 0. This is important since c is not necessarily continuous at 

X- = 0. 

iHere we used P'"' J dx-A'l{x~)d+A'r{x- ) = i/°^'=7f7f + rfa;-i^(x-)a+i^(x-) which foUows from 

eq-dn?!). 
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3.2 Solving the constraints 



Our strategy now is the fohowing. We should solve the two constraint equations, eqs. (|3.16| , p.l8| ) 
and determine the commutation relations of the non dynamical field 7. Then we must substitute it 
back into the Hamiltonian and express the Hamiltonian in terms of the canonical degrees of freedom 
A. 

We will do so by expanding the constraint equations and the Hamiltonian in powers of g. When 
doing so we must have some knowledge of the parametric dependence of the valence charge density 
j on the coupling constant g. The expansion in principle can be performed for any parametric 
dependence. In this section we take j to be of order as in the JIMWLK limit, and will collect 
all contributions to the Hamiltonian of order 1/g and order 1. It turns out however that this same 
resummation collects the leading terms in g also for any g < j < 1/g- We will discuss this point in 
detail in the discussion section. Thus even though in this section we treat explicitly j as being of 
order 1/g this should not be construed as limiting our calculation to the JIMWLK limit. 

Thus our aim in this section is to expand 7 to 0(1), obtain the Hamiltonian to 0(1) and 
diagonalize this 0(1) Hamiltonian exactly. Further corrections to this calculation are strictly per- 
turbative (small corrections in powers of g for any parametric dependence of j on g) and will not 
be considered here. 

To order 1/g the operator 7 satisfies the 'classical equations' 7" = 6": 



di m{x) - ht{x) - gr'' h\{x) m{x) = . 



(3.21) 



To this order the commutation relations are calculated as 



5f{z) 



-ig 



D, 



dD 



5f{z) 

{x,z) r'^'fiz) 



db 



(3.22) 



where D is the transverse covariant derivative in the 'classical' background field b: Df^ = di5'^^ — 
gjacbjjc^ Eg ( ^|3 22| ) is the leading order result in g. Note however that it is exact in the weak field 
limit, where the field b is linear in the valence charge density j. 
Eq. (|3.22|) can be further simplified, using the identity 



gr'^'fiz) = gr'^d.btiz) 

Thus finally, to leading order in g 



[d^,d^-D;\"' = {dD - Doy'. 



(3.23) 



bt,^] = [bl A)]=0. 



d'^{x,y)^i[D,—D, - a,^D,rix,y) 



ab I 



dD 



Dd 



(3.24) 
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Note that although 7 itself is of order 1/ g, the commutator of two 7's is of order one. It is 
thus clear that we will not need higher order corrections to the commutator eg. ( ^.241 ) in the 0(1) 
calculation. 

To order 0(1) we write 



K + Q 



(3.25) 



where C satisfies the equations: 



ab/-b r^ab/-b g 



0), 



The solution to these two equations is easily found as 

1 



D. — id-D)Aix- 



0) 



(3.26) 



(3.27) 



where the product on the right hand side is understood in the matrix sense over all indexes (including 
transverse coordinates). Note that the ordering of different factors of b in eq.( |3.27| ) is irrelevant, 
since the covariant derivative involves gb, and the commutator of two such factors is 0{g'^) and is 
thus of higher order than the one we need to keep. 

The canonical structure to 0(1) follows from Eqs.( p.25| , p.27D 



bt{x),A]{y)] = [Q{x),A]{y)] = -ie{y 
[c-(x),7,'(y)]=^6(x-) 



ab 



ab 



[ct{x),c%y)] = --e{x--y-)6t,'{x-y) - -e{x-)e{y-)C^^{x,y) 



where for future convenience we have defined 



Oj(a:,2/) = |29,-i^D,-2A^9,j 



ab 



{x,y) . 



(3.28) 



(3.29) 



3.3 The Hamiltonian and the equations of motion to 0(1) 

Next we express the Hamiltonian to 0(1) in terms of the field c". 
For the magnetic piece to 0(1) we have: 

Gt,{e{x-h + c) = e{-x~)[d,c, - d,c,] + e{x')[D,c, - D,c,] + 0{g) . 

For the electric piece, using the constraint and after some algebra, we obtain: 



Il- = -^[D.d^A-j6ix-)r 



[e{-x-)d,c, + eix-)D,CiY + Oig). 



(3.30) 



(3.31) 
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All said and done the Hamiltonian to 0(1) is 
H = - ^ j dx-(fx [e{-x-) c1{x~,x) c^(x-,x) + e{x')c1{x-,x) D^""^ c\{x-,x)\ . (3.32) 

This is the Hamiltonian that we have to diagonalize. The most efficient way of doing this is first 
to find the complete set of solutions of classical equations of motion, and then expand the quantum 
field operators in the canonical creation and annihilation operators with the coefficients given by 
the solutions of classical equations. The classical solutions have to be properly normalized in order 
that the quantum field operators satisfy correct commutation relations. 

We start by deriving the equations of motion. Using the commutation relations eq.( p.28|) we 
obtain 

id^d-ct{x) = [H,d+cnx)] = I dy-[d+ct{x),c%y)] [e{-y-)d' + e{y-)D']';'^cl{y) 

= -l[0(-x-)9V + e{x-)D^''Yi{x) 

-'-6{x-)C:;'{x,y) I dy- [-e{-y-)d' + e{y-)D']';;^cl{y- ,y) , 

(3.33) 



where C°'j' is defined in ( ^3.291 ). Integrating these equations (avoiding the singularity at y = 0) 
gives 



0- ro- 



dy-d'c{y) = i / dy-d+d-c{y) = ia"c(0") , (3.34) 

) J — oo 

• poo poo 

- / dy-D^c{y) = i \ dy-d^d'ciy) = -i9"c(0+) . (3.35) 
2 Jo+ Jo+ 



10+ Jo+ 
The last term in eq.( p.33|) can be rewritten as 



-'-6{x-)Q^{x,y) J dy- [-e{-y-)d' + e{y-)D']'^cl{y) = -t6{x-)C^;>{x,y)d-c'^{0) , (3.36) 
so that finally the equations of motion are 

t[d+ + 6{x-)C]^;{x,y)d-c%y) = -'-[e{-x-)d' + e{x')DY{x,y)c'i{y) . (3.37) 
Matching the discontinuity across x~ = gives the relation 

c^(0+,x) -c^(0-,x) = -l.CtHx,y)[c^{0\y) + 4{0-,y)]. (3.38) 
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The solution to the equations of motion can be written down explicitly. At negative x~ this is just 
a free equation, and thus the solution is a superposition of plane waves. At positive x~ the solution 
is again a superposition of gauge rotated plane waves. This can be written as 



exp{zp x'^} j d^q 



_ ■_ _ _ 

e{-x )exp{i—x }vl^g{x) + e{x )exp{i—x jv^t^ix) 



3.39) 

Except at x" = this solves the equations of motion with given p~ for arbitrary f Here q is the 
degeneracy index. In the free theory the index q would stand collectively for transverse momentum 
k, polarization index i and color "polarization index" a. In the present case q also stands for i and a 
as well as some continuous degeneracy. For simplicity of notation we will not differentiate between 
discrete and continuous parts of q. In the following, integral over q stands both for the integral over 
continuous part with appropriate measure as well as for summation over the rotational and color 
'polarizations'. 

Eq. dpI) imposes the condition 



vrix)-vrix) = -lct;'ix,y)[v',^y)+vtiy)] . 
This equation can be equivalently rewritten as 

where the projectors T, L, t, I are defined as 



T ab 



ab 



ab. 



li 



6ij lij 



Eq.( p.41|) is solved by 

for arbitrary Vj. Thus we can write the solution eq. (|3.39| ) in terms of one set of functions 



p 1^ 



(3.40) 
(3.41) 

(3.42) 

(3.43) 

fx) as 



exp{ip x^} I d^q 



52 

e(-x )exp{«— X ][t-l]ijvl.^{x) 



+e(x )exp{i— X }[T-L]ijV_^(x) 



(3.44) 



On the level of the classical solution, the normalization of the functions ^^^^(x) is arbitrary. How- 



ever, in order to use eq. (|3.44|) as the basis for expansion of the operators c in terms of canonical 
creation and annihilation operators the normalization of ?^pi^(x) has to be determined. This will be 
done in the following subsection. 

As a corollary to this subsection we note that the classical field h does not commute with the 
Hamiltonian and is therefore not constant in time. Calculating the commutator we obtain 



id- mix) 



dy~mx),c%y)] [e{-y~)d' + e{y-)D% cl{y) 



i4?(y,x) J dr [-e{-y-)d' + e{y-)D^]';^ cliv). 



(3.45) 
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Using eqs.( [3.34P-35|) this can be written as 

zd~btix) = 2zd~ \^D,^D, - A^/^,j 4(0) . 



(3.46) 



This can be interpreted in the following way. Let us define the operator b, so that it has the 
same exact matrix elements on the valence part of the Hilbert space as b, but commutes with the 
operators c. Then we can write 



1 



bt = bt + 2\D,^D,-a,—D,\ c^iO), 



dD 



ab 



(3.47) 



and 



li 



ab 



(3.48) 



This form will be convenient for calculating correlators of 7 in the vacuum state. 
3.4 Normalization of the eigenfunctions and the vacuum state. 

Given that the 0(1) Hamiltonian is quadratic, and having found the complete set of solutions of 
the classical equations of motion, we can find the quantum vacuum state. 

The vacuum state of the Hamiltonian eq. ( |3.32| ) is the Fock vacuum of the canonical annihilation 
operators Pp-,q defined in terms of c by 



dp 



+e(x-)e^i^^"[T-L]^?(x,i/)t;Jl^(y) 



(3.49) 



where the integral over the transverse coordinate y is understood but not written explicitly. The 
operators (3 satisfy canonical commutation relations 



= {2 7r)6{p--p'-)6{q-q') 



(3.50) 



Existence of such a set of canonical operators is guaranteed if the set of solutions of the classical 
equation is complete and the functions v entering eq.( p.49|) are properly normalized. To find the 
correct normalization of these functions we require that c satisfy eq.( 3.28| ). 

We concentrate on negative x~ and y~ first, so that only the first term in the sum in eq. ( p.49| ) 
is important. For simplicity we suppress the color indexes and also the factor t — I, thus we are 
working in terms of v~ rather than v. Consider the commutator 



[c'ix),c^iy)] 



dp 


r - 
li — 

e 2p" 


S2 

°y 

2p- 


271 






dp~ 


' r - 
li — ^x 

e 2p- 


32 
y 

2p- 


27r 







y^} 



V 1 *_"' iv) 



— 2 [ 3; ^ u \ 

g 2p- 2p- " ' 



v*J (x)v f iv) 



(3.51) 
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We have defined the 'correlator matrix' 



(3.52) 



Note that this matrix fully determines the commutators of c, and there is no need to find the 
individual functions fp- g. Different choices of the functions v which give the same W correspond 
to unitary rotations of the set of the canonical operators f3. 

To determine the correct normalization we first note that taking 



z p 



(3.53) 



would give canonical commutation relations for the fields c. With this expression for W we can 
change variables p~ —p~ in the second term of eq. (|3.51|) to get 



dp 



4:71 P 



-e 



2p- 



-y }x2 



6\x - y)6'^ = — 6'^6\x - y)e{x- - y-) 



(3.54) 



where the last line follows by change of variables p —>■ d'^/2p . To get the e-function in the 
commutator we have to regulate the singularity in in the symmetric way 



1 

p- 



1 

p- 



+ ie 



+ 



(3.55) 



To reproduce the extra term in the commutator of c* (the second term in the last line of eq. (|3.28|) ) 
we modify the matrix W in the following way 



_ ^) + ^C*^(x,y)] + ^^\b^^\x -y)- \c'\x,y)\ 

— + I — — %t I 

p- p- 



(3.56) 

The new term we have added is imaginary and even with respect to p~ — > —p~ . Thus it is still true 
that the two terms in eq. (|3.51|) are equal. The extra term under the change of variables p~ — >■ l/p~ 
gives 

and thus generates the term in the commutator independent of x~ and y~ . The result is precisely 
the last term of eq.(F^). It is a matter of some straightforward algebra to check that with W 
defined in eq.( p.56| ) the correct commutator of the fields c is reproduced also for other values of x~ 
and y~. The following identities come handy in this calculation 



-C 



(T-L)(t-/); 



1--C 
2 



{t-l) 



^ + 2^ 



(T-L)- 



{t - l)C{t - I) 



-c- 



{T - L)C{T - L) = -C. 



(3.58) 
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Returning from v to v we conclude that the operators (3, f3^ in the representation eg. ( p.49| ) 
have canonical commutation relations when (we use eg. (pi 



— + ie — ie 

p p- 



p p 



5'^^{x - y) 



(3.59) 



We thus conclude that the vacuum of the Hamiltonian eg.( p.l|) to 0(1) is the Fock vacuum of 
the annihilation operators (3 related to the original gluon field operators through 



A1{x ,x) = e{x ) 



b^(x)-2D,^{d-D)ix,y)c{0,y) 



+ ctix ,x) 



eix 



1 



1 



btix) + 2 \ D,—D, - D,—dj \ (x, yy^iO, y) 



dD 



ab 



+ c^(x-,x) (3.60) 



with the field c^(a; , x) expressed in term of (3 and in eg. (|3.49|) with the normalization eg. (|3.59|) . 
This completes the diagonalization of the light cone Hamiltonian to 0(1). 



4. The Bogoliubov operator 



The calculation of the previous section can be viewed as the diagonalization of the light cone 
Hamiltonian. Although we have only found the vacuum state, guite generally the diagonalization 
is affected by the action of some unitary operator Q. Namely for the case of a guadratic operator 
H 



(4.1) 



P~,Q 



The explicit knowledge of the operator Q, or alternatively the knowledge of its action on all the 
degrees of freedom of the theory furnishes much more information than just the vacuum wave 
function, as it also in principle can give us the wave functions of excited states, which are necessary 
to calculate more exclusive properties than the forward scattering amplitude. The aim of this 
section is to find explicitly the action of Q on the degrees of freedom of the theory. 

Part of the answer to this guestion is already furnished by eg.( |3.6CI|) which can be viewed as 
the transformation of the vector potential if we read the left hand side as Q^AQ and the canonical 
operators P and P^ in c on the right hand side as the original gluon creation and annihilation 
operators a and a"^. The missing piece of information is the transformation of the valence charge 
density. This is the guestion we address now. 

First, it is clear from eg. (|3.60|) that the transformation is of the Bogoliubov form, namely 



Q = C B = exp 









EA 


exp 


-AM A 
2 



(4.2) 
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where E and M are operators which depend on the charge density j but do not depend on the 
soft fields A. We do not indicate exphcitly the indexes and coordinate dependences of E and M 
for simphcity. Those should be clear from the context. Here C is a purely coherent state operator 
- exponent of an operator linear in A, while B has no linear term in the exponent. The coherent 
operator is easy to find by inspection, since it is the only one that induces the shift of the soft field 
(the very first term in eq. (|3.6CI|) ): 



C = exp 



2i / d^xb'^{x)A1{x- = 0,x) 



(4.3) 



The Bogoliubov part of the transformation, the operator B is more difficult to determine. Rather 
than looking for the explicit form of the operator B in terms of j, we will find its action on the 
degrees of freedom of the theory by considering sequential action of C and B on A and matching it 



onto eq.( |3.6CI| ). 

It is important to remember that we need to know the transformation of the color charge density 
only to 0{g). Only this order contributes to the JIMWLK evolution as explained in detail in [p!0| . 
Thus we will determine the action of B on the fields to this order only. 

We first note the following 'combinatorial' identity. For any operators O and L 



e-^Oe^ = 0+[0,L] + ^[[O, L], L] + i[[[0, L], L], L] + .. 



(4.4) 



Using eq. 



we have for C of eq. (|4.3|) 



C^At{x)C = At{x) + e{x-)b^{x) + e(x-) dt';{x,y) A%y- = 0,y) 



+ ^e(x-) / [d^^{x,y),bl{z)]A%y- = 0,y)Al{z- = 0,z) 

>J Jy,z 

C^r{x)C = r{x) + 2^|(^9Z}-i^-ljD,j {x,y)A%y- = 0,y) 



+ 2i 



ab 



{x,y),bl{z) 



(4.5) 



A%y- =0,y)Al{z- = 0,z).iA.6) 



To find the action of the Bogoliubov operator, we imagine diagonalizing the Hamiltonian first 



by acting with C and then subsequently acting with B. Transforming the Hamiltonian eq. (|3.1| ) with 
C obviously leads to 

H[A,j]C ^ H'[A,j] = HpACC^jC]. (4.7) 

It is straightforward to see using the expression for the transformed fields eqs.( |4.5| , [4.6[ ), that if we 
substitute for A in the function H' the following expression 



i^(x)^c^(x) + eix')A^';ix,y)c%y~ = 0,y) 



(4.8) 
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with 

r 1 1 1 

Ag'(x, y) = D,—D, + D,—D, - 2D,—d, (x, y) , (4.9) 



we obtain to 0(1) precisely eq.( 3.32 ). This substitution should be equivalent to the action of the 
Bogoliubov operator 

H'[A, ]] B = H'p AB, B^ j B] = H"[A, j] . (4.10) 

In other words, up to (and including) 0{g) terms the action of the Bogoliubov operator B on the 
field A is 

At{x,j) ^ B^A1{x)B = ctix) + eix-)A1^ix,y)c%y-=0,y), (4.11) 

where the field c on the RHS is understood as expressed in terms of the canonical creation and 
annihilation operators a and (rather than f3 and /?^)^. 

Our aim is now to find the transformation of the color charge density under the Bogoliubov 
transformation which induces eq. (|4.11 ). This is indeed possible, even though we do not know the 



explicit form of the operator B itself in terms of the fundamental fields. The key is given by the 
following chain of arguments. 

Consider a general Bogoliubov operator of the form 



B = exp 



-Ai Mij Aj 



(4.12) 



Here we denote all indexes/coordinates of the field A by a single index i. The fields A are assumed 
to satisfy the commutation relation 

[A, Aj] = P,, (4.13) 

with some matrix P. Quite generally the matrix M is symmetric and anti hermitian, while P 
is antisymmetric. The matrix M depends on the charge density and the coupling constant only 
through the combination g j. 
Consider the transformation 

Ak = B^AkB = Ak + iPMA)k + ^iPMPMA)k + ^iPMPMPMA)k + ...= [e''%iAi. (4.14) 



Here we have used the identity eq. ([4.41). Also, consistently with our counting of powers of the 
coupling constant we have neglected all and any terms involving commutators of gj which enter 
into M, since each such commutator brings a power g"^. 
Now to order 0{g) we have 

[f,M,,]=tgr''r^. (4.15) 



2 



We note that stric tly s peaking to make this identification we should also substitute into H' the transformed 



expression for j in eq.(4.10), which we do not know at this point. However as we will see below and is simple to 



understand by straightforward counting of powers of g, the operator B induces transformation of j only to or der g 



Since we only need the Hamiltonian to 0(1) it is therefore perfectly consistent to keep j unchanged in H' eq.( 4.1C ) 
for the purpose of the identification of the Bogoliubov transformation of A. 
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Thus consider the transformation of j"'{x) induced by the action of B in eq. (|4.12 ): 



(4.16) 



+ i^PMPM + MPMP^ - 2MP^PM ] A + ... 
3! \ dj" oj" dj" 



Here again we neglected all commutators of gj in M beyond the first term, as they are all higher 
order in g. The negative signs come from transposing the antisymmetric matrix P. We can now 
check explicitly that eg. ( [4.16|) is expansion in powers of M of the following expression 



r + ^gr'rAP-'^j^^ 



(4.17) 



with A defined in eq.( ^4.14 ). Remembering that in our case P = |e(x —y ) whose inverse is —2id^, 
and restoring all the indexes and coordinate dependences we obtain 



5f{x) 



(4.18) 



with A given in eq. (|4.11 



An equivalent way of obtaining this result is to require that the transformed fields satisfy the 
same commutation relations as the non transformed ones, the transformation being unitary. Using 
the explicitly known commutator of the field c one can easily show that 



[A^{x-,x),A'{y~,y)] = --eix' - y-)6t ix - y) 



(4.19) 



[rix),A'iy-,y)]= l[rix),fiz)] 



tgrYix) 



SA'^j{y-,y) 

6f{x) 



In this expression we should understand ^ as a function of j at fixed a. It is easy to check that 
with the transformation eq.( [4.18|) to 0{g) 



[rix),A%y-,y)] = 0. 



(4.20) 



This is straightforward after noticing that the last term in eq.( [4.18| ) can be written as 



gr'fix) / dy-d'yd+A''{y-,y) 



SA';{y~,y) 

6f{x) 



dy-d'yd'zA''{y-,y)[f{x),f{z)] 



sA'^{r,y) 

5f{z) 



(4.21) 

Therefore we conclude that the transformation eqs.( |4.11|j4.18| ) does indeed preserve canonical com- 
mutation relations of the fields. 
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We can now put all the elements together and write down the transformation that the operator 
fl induces on the fields: 



^l^'AKx', x)n = Ci{x~, x) + e{x~ 



ht{x) + 2 I \ D,—Dj - D,—d, \ (x, y)cp, y) 



dD 



ab 



-ex 



2i, 







/f 





ab 



ixz)n\z) I dy-d+A^iiy-,y) 



SA''^{y-,y) 



+j[dt;{x,y)Mz)]c%0,y)cl{0,z) 



n^rix)n = fix) + 2^1 (^dD^ - 1^ ix,y)c%0,y) 
+9rY{x) ! dy-d'yd^A'Ay-,y)- 



(4.22) 



5f{x) 



+2i 



y,z 



1 



dD—-l]D^} {x,y),bl{^) 



c'My)cl{0,z). 



(4.23) 



Here A is given by eq. (|4.11|) and the field c is understood as expressed in terms the canonical 
creation and annihilation operators (3 and as in eq.( |3.49D . The first line of eq.( |4.22| ) coincides with 
eq. (|3.60|) . The second and third lines are the 0{g) terms. They are given here for completeness even 
though they do not contribute in the calculation of the previous section and also do not contribute 
to the transformation of the soft color charge density eq. ([4.25|) . 

Eqs.( [4.2^ , [4.23[ ) are the main result of this section. They give the explicit action of the diago- 
nalizing operator Q on the fundamental degrees of freedom of the theory. 

Finally, for completeness we give the expression for the transformation of the total charge 
density. This is the observable directly relevant for the calculation of the scattering amplitude. It 
includes the contribution of the valence and the soft modes 



r{x)=rix)+gr'^ I dx-Atix)d+A^ix) 



(4.24) 



Collecting the formulae given above we find 

fitj" (x)^] = f{x) + 6ir{x) + 62r{x) 



(4.25) 
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with 



Sif(x) = 2 



SDj^^D, - a, 



ab 



{x,y)c%0,y) 



+9r'f{x) j dy-d^A%y-,y)^^, 



with A given by eq.( [4.11|) . Here 



Sb'^ix) 



dD 



ab 



{x,y) 



(4.26) 



(4.27) 



(4.28) 



As a consistency check with the calculation of the previous section we note that eq. ( |4.26| ) coincides 
with the divergence of eq.( |3.48| ). 



5. Reproducing JIMWLK/KLWMIJ. 

As a cross check on our derivation we reproduce in this section the two known limits of the high 
energy evolution - the JIMWLK evolution equation (the high density limit) and the KLWMIJ 
evolution equation (the low density limit) . 

5.1 The JIMWLK kernel 



Under boost the color charge density j transforms into J of eq.( |4.25| ). To derive the evolution of 
the functional W we have to calculate the correlation functions of J over the soft gluon vacuum, 
that is over the Fock vacuum of operators f3. In the JIMWLK limit it is only necessary to know 
two correlators, 



■x"-^{x,y) = lim 



(0/^1 h3''{x)5i3\y) |0^) 



a''(x) 



lim 



(0^1 ^j-(x) |0^) 



(5.1) 



since Sij ~ gj and 823 ~ and so only these two correlators contribute to the evolution of 
{j{xi)...j{xn)) to relative order g"^. In fact our task is somewhat easier, since we can avoid the 
calculation of {Sj2{x)) using the following argument. In terms of x and a the evolution kernel has 
the form 

-^^^^"^ = y\x,y) ^i^, , + a"(x)-A- . (5.2) 



63'^{x) 63\y) 



53^{x) ■ 



However it was proved in that the evolution kernel has to be a Hermitian operator (on the 
space of functions of 3). In conjunction with the fact that (j"'{x) is real, since it is a diagonal matrix 
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element of an Hermitian operator (on the QCD Hilbert space), it means that a is rigidly related to 
X so that the evolution kernel is 



E 



JIMWLK 



1 5 



2 53-{x] 



(5.3) 



This property of the JIMWLK kernel is of course well known and has been first noted by Weigert 
in the last reference in [|1^]. Thus our task is first to calculate x^'^i^)!!) ^"^^ then to show that 
the resulting evolution equation is equivalent to the standard form of JIMWLK which involves 
derivatives with respect to the unitary matrices U rather than with respect to the charge density j. 
We start with the calculation of Xi defined as eq.( |5.1|) . In preparation we calculate 



(0^|c^(0,x)c;(0,y)|0;3) 



8Ar 
1 



^^1 [t _ / + T - L] [t - / + T - L] J (x , y ) 



2'Kp 



An 



[l-l- L + IL + Ll]t'{x,y). 



(5.4) 



Using eq.( [4.26|) we then find 

(5ij-(x)5ij^(y)) 4 



Ar 



1 



dD\ 



1 



Ar 
1 1 



dD—Di - 
Do 



{x,u){ct{Q,uyMv)) 



d.-D.—Dd 



db 



92 



ab 



{x, y) . 



{v,y) 

(5.5) 



5.2 From j to U. 



To get the evolution equation in the familiar JIMWLK form we need to change variables from j to 
the single gluon scattering matrix U. The matrix U is defined as the matrix of the two dimensional 



gauge transformation which transforms the 'classical field' b to zero value |jTO 

ab 



U-\x) = !^Vexp[zg jjy.T'^b'riy)] 



(5.6) 



where the contour C starts at some fixed point at infinity in the transverse plane and ends at the 
point X. The matrix U does not depend on the curve C but only on its end point, since the field b 
is two dimensionally a pure gauge. Using this definition we have 



9 / dyi 



c 



U{x)U\y)T'fMu{y) 



-1 ab 



g j dy,[U{x)U^{y)T'^U{y)]'''[D,^ 
[U\y)T''U{y)t = T^,U'^%) . 



Now we use the identity 



dc 



(5.7) 
(5.8) 
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Substituting this into eg. ( |5.7|) , and using the fact that J^dyidiF{y) = F{x) we find 



UT 



dD 



{x,z) . 



(5.9) 



This makes it possible to rewrite the real part of the JIMWLK kernel in the following form 



6 6 



x,y 



6j-{x)6f{y) TT J,^y6U-^{x)6U<^^{y) 



6 



6 



[U{x)TY[Uiy)T 



diem 



1 1 Ln/)1 

g2 + JJ2 Q2^ ]J2 ^ Q2 



Im 



Now remember that 



[U{x)T 



bial 



-Tr 



6W{x 



-U{x)T^ 



-J 



R 1 



(5.10) 



(5.11) 



where Jr is the operator of right rotation on matrix U. We also note that 



rj 1 / \ 1 'f-j yi 



1 



1 



D, — {x,y) = —U\x) 



Now, using eg. ([5.3[) we can write the complete kernel as 

JIMWLK f — z)i{y — z) 



{x - yf 



u{y)- 



(5.12) 



H 



27r2 



x,y,z 



(x — zy{y — zy 



[Jl{x)Jl{y) + J^(x) J^(y) - 2Jl{x)U^\z)j'M] (5.13) 



with J£(x) = f/"''(x)J^(x). This is by now one of the standard forms of the JIMWLK kernel, 
see 1311. 



5.3 The KLWMIJ evolution 

Although our derivation has been formally in the high density limit, as we noted in the introduction 
and as we explain in the next section the result egs. (|4.25|J4.2^ , [4.27] ) is in fact valid for all physically 
interesting situations, including the low density case j = 0{g). For the low density case we have 
to reproduce the KLWMIJ evolution eguation It is easy to see that this is indeed the 

case. Examining the action of the Bogoliubov operator B on the fields, we see that in the weak 
field limit they are sub leading. The shift of the vector potential affected by the coherent part of 
the operator C is of order b j, while any correction introduced by B is of order gb ~ gj. This 
is also true in the strong field case, however for j ~ 1/g the corrections due to B are 0(1) and 
therefore could not be neglected. In the weak field case these are not only sub leading but also 
genuinely perturbative! We can therefore neglect the action of B altogether. Thus in this limit the 
operator Q reduces to the coherent operator C with the 'classical field' b given by the leading order 
perturbative expression. This is precisely the operator that was used in |2^ to derive KLWMIJ 
evolution eguation. Obviously, repeating the same derivation we obtain the same result. 
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One important thing to be noted here is, that in order to derive KLWMIJ we are not allowed to 
expand the correlators of the transformed charge density to first order in 62 j eq. (|4.27|) as is done to 



derive JIMWLK equation. The reason is very simple. When j ~ 0{g), the second term on the RHS 
of eq. (|4.27|) is of the same order as j itself . Therefore its contribution to the evolved correlators of 



J has to be resumed to all orders. This is indeed what is done in the derivation of [^. It is the 
resummation to all orders in f°'^'^J dx~A''[x)d~^A'^{x) that is responsible for the appearance of the 
'dual Wilson line factor' 

R{x) = exp T" 



6j^{x) 



in the KLWMIJ evolution equation [24]. 
6. Discussion 

In this paper we have carried through the diagonalization of the QCD light cone Hamiltonian in the 
presence of a valence charge density j. We found that for large valence charge density, to 0(1) the 
vacuum is the Bogoliubov transform of the free gluon vacuum. We have also found the action of the 
Bogoliubov operator on the dynamical variables of QCD including the valence color charge density. 
The evolution of hadronic wave function to high energy increases the longitudinal momentum of 
the gluons in this state. Thus more energetic gluons scatter on the target leading to the evolution 
of the hadronic scattering matrix. 

6.1 The calculation is valid for any j. 

We have shown that when the valence charge density is large j ~ 0{l/g) the wave function we 
found leads to the JIMWLK evolution equation. However, our calculation itself is valid beyond 
the high density limit, and does in fact give the leading solution of the light cone Hamiltonian for 
all physically interesting magnitudes of the color charge density j ~ 0{g'^)\ — 1 < n < 1. The 
precise statement is that relative corrections to the solution we have given here are proportional to 
a positive power of g at any interesting value of the valence charge density. 

To see this, recall that the basis of our approach was the perturbative solution of eqs. (|3.16|J3.19|) . 



We have solved eq.( 3.1(j| ) exactly, while eq. ( |3.19 ) was solved treating the second, third and fourth 



terms on the left hand side (LHS) as perturbations. The solution of this pair of equations to leading 
order in the coupling constant is always of order 7 = = 0{g^) for — 1 < n < 1 . The magnitude 
of corrections is easy to estimate. Since, by definition the field c is 0(1), we have 



^r^{7f(x),c^(x,0)} = O((7^+"), 

gr"^ [ dx'4{x')d~'c^{x~) = 0{g) . (6.1) 



The first term is always smaller than the second. It always scales as a positive power of g and 
therefore can always be treated perturbatively. 
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The second term is also small as long as n 7^ — 1. It is a factor g smaller than the zeroth order 
solution and thus again can be safely treated perturbatively. For n 7^ — 1 it can be neglected since 
its magnitude is a positive power of g. The case n = —1 is a bit different, since then this term is 
0(1) and so has to be taken into account, which is what we did above. 

Finally the third term is always 0{g). It can be neglected for all n 7^ 1. For n = 1 this term 
is of the same magnitude as j and thus it may seem that it has to be taken into account already 
in the leading order. However this is not the case for the following reason. The vacuum of the 
Hamiltonian of the field c at j ~ 0{g) is a free vacuum. This state is annihilated by the 'soft' 
color charge density operator f"'^'^ j dx~ (^{x)d^ c^{x). Thus this operator only gives non vanishing 
contribution to 7 in the sub leading order in g, where the vacuum is not a free vacuum anymore. 
Thus we see that for all —1 < n < 1 our solution of eqs.( p.l€ , 3.19 ) keeps the leading terms and 



for n = —1 also the important sub leading term of 0(1). The terms that we omit are not only 
suppressed by a positive power of g relative to the terms we keep, but also vanish in the limit g ^ 
at any j. 

To reiterate, our procedure keeps all the terms that are important for physically interesting 
values of the color charge density. This is not to say that our solution can be considered as a 
leading order of some expansion which has the same expansion parameter for all n. The corrections 
to the leading term may have different magnitude for different values of n, and thus the properties 
of the expansion are different at different values of n. At this point however we are not interested 
in the sub dominant corrections and will not discuss this issue any further. 

6.2 What JIMWLK misses? 

We want now to return to the point briefly mentioned at the end of the previous section. Even 
though our diagonalization procedure and the solution for the vacuum wave function is valid for 
any j, the derivation of the evolution equation for the scattering amplitude involves one extra step, 
and that is adding the charge density of the soft gluons to the valence charge density. For n 7^ 1 
this is a perturbative proposition, since the soft gluon charge density is parametrically smaller than 
j itself. Thus for the derivation of the JIMWLK evolution equation one expands to first order in 
the soft gluon color charge density, the second term on the RHS of eq.( [4.27| ). For the KLWMIJ 
evolution on the other hand all powers of the soft gluon color charge density are resummed. The 
addition of the soft gluon charge density is achieved by acting on any observable function of j by 
the shift operator of the form 



Ra = exp 



(6.2) 



In the KLWMIJ limit only one gluon is produced at one step of the evolution with probability 
of order a^, and thus jgofti^) = when acting on the component of the wave function which 
contains this extra gluon. The phase factor of eq. ( |6.2|) therefore simply becomes the dual Wilson 
loop R. In the general case however the action of the Bogoliubov operator B produces an arbitrary 
number of gluons. For j ~ 1/g the number of gluons of order 0(1) is produced with probability of 
0(1), while with probability 0{g) one can produce 0(1/ g) extra gluons. The phase factor becomes 
a product of dual Wilson loops R{xi)...R{xn) when acting on a component of the wave function 
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with n extra gluons. Now the JIMWLK equation is vahd when a large dense target scatters off a 
small perturbative projectile. In this situation each gluon in the target wave function undergoes 
only a small number of scatterings on the projectile. In fact the leading order scattering on a 
small target is only due to two gluon exchange. This corresponds to expansion of each dual Wilson 
loop factor R to second order in S/6j. It is also true that in this situation it is unlikely that two 
or more produced gluons scatter simultaneously. Indeed the expansion of eq.( |6.2|) is equivalent to 
approximating the scattering amplitude of the configuration of n produced gluons by the sum of 
the individual scattering amplitudes. 

Recall that in calculating the evolution of any correlation function of j in the JIMWLK ap- 
proximation, we only keep terms of the first order in 62 j and of the second order in 61 j. Thus the 
correction to the scattering matrix S = exp{i j ax} due to the evolution is at most of second order 
in the target field ar- This is another way of saying that the whole system of soft gluons produced 
in one step of the evolution scatters on the target only via the two gluon exchange. To be a little 
more precise we have to remember that while calculating the evolution of the scattering matrix, a 
factor R accompanies not only each soft emitted gluon but also every factor of j in the operator Q. 
The eikonal scattering matrix of the projectile wave function on the target field ay is given by (we 
drop the transverse coordinate dependence to simplify the notations) 

= (vl/[j1|fit[j,a,at]e*(^-"+^-"°/>-fi[j,a,at]|vl/[j]) = (M/[J]|^]t[^■, a]e^^-"-fi[j, i?a, i?at]|M/[j]) 

= mj]\n^[j,a]n[Rj,Ra,Ra^y^'''^"-mj]) , (6.3) 

where l^&fj]) is the valence wave function and the functional derivatives in R act only on the eikonal 
factor e*-'""!'. The first equality is the reflection of the fact that multiplying every soft gluon 
creation operator by R is equivalent to shifting the charge density j by the charge density of this 
soft gluon. The second equality follows from commuting of the operators j in Q with those in the 



eikonal factor as explained in detail in ||3^. Since every R becomes an eikonal factor after acting 
on e*-'""!', multiplication of j by R in the second line in eq.( |6.3| ) physically corresponds to the effect 
of scattering of the valence charges involved in the emission of soft gluons. Thus the expansion of 
all the factors of R to second order in jj approximates the interaction of the whole system of soft 
gluons emitted in one step of the evolution plus the valence charges involved in their emission (in 
the following we will refer to this system as "soft gluons" to avoid lengthy and wordy descriptions), 
with the target by a two gluon exchange. 

The JIMWLK evolution therefore does not take into account multiple scattering corrections to 
the amplitude due to simultaneous scattering of two soft gluons emitted in the same step of the 
evolution. This is not to say that the JIMWLK evolution does not allow any multiple scattering 
corrections at all. In particular the probability that a soft gluon scatters simultaneously with some 
of the valence gluons not participating in its emission, is accounted for. We will refer to these 
multiple scattering events as "long range multiple scatterings" to emphasize the fact that the two 
objects that scatter simultaneously have vastly different rapidities. This as opposed to "short range 
multiple scatterings" where both objects have similar rapidity, which are taken into account by the 
KLWMIJ evolution. 

It is these long range multiple scattering corrections that unitarize the scattering amplitude 
in the JIMWLK approximation. If no multiple scattering corrections where included at all, the 
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amplitude would not unitarize even though the coherent effects in the wave function are taken into 
account exactly. Recall that the charge density itself does not saturate even in the dense regime, 
although its growth with rapidity is much slower than in the BFKL approximation ||2^ . In particular 
in the BFKL (or equivalently KLWMIJ) limit the color charge density grows exponentially with 
rapidity 

f(F)(xj(0)V^, (6.4) 



while in the " saturated regime" the growth is a random walk process and thus 

f ocj{Of + kY . (6.5) 

Since the charge density does not stop growing even in the saturated regime, the scattering am- 
plitude would not saturate if no multiple scattering corrections are taken into account. It is thus 
precisely the long range multiple scattering corrections that stop the scattering amplitude from 
growing beyond one in the JIMWLK approximation. 

Eqs.( |6.4D and (|6.5|) in fact clearly indicate that the short range multiple scatterings are dominant 
in the KLWMIJ regime while the long range multiple scatterings are dominant in the JIMWLK 
regime. Consider first the evolution of a dilute projectile (KLWMIJ evolution). According to 
eq.(|6.4|) the color charge density grows exponentially fast and is always (at large enough rapidity) 
dominated by gluons created in the last rapidity interval of the size AY ^ -. Thus the dominant 
multiple scattering effects indeed are due to the simultaneous scattering of two or more gluons at 
approximately the same rapidity - the " short range multiple scatterings" . On the other hand in the 
JIMWLK regime where eq.( |6.5|) is valid, the color charge density is uniformly distributed in rapidity. 
Thus clearly the dominant multiple scattering corrections are due to simultaneous scatterings of 
gluons at far away rapidities - the "'long range multiple scatterings"'. 

We thus see explicitly that while the KLWMIJ evolution takes into account all multiple scatter- 
ing effects but does not include nonlinearities in the evolution of the wave function, the JIMWLK 
evolution fails to take account of the short range multiple scattering corrections to the amplitude. 

6.3 Short range multiple scattering and the dipole - dipole amplitude. 

In relation to the preceding discussion we want to comment briefly on one aspect of the Pomeron 
loop correction to the JIMWLK evolution. In particular recently much attention has been devoted 
to scattering of two unequal size dipoles. In this context there has been much discussion of the 
effects of discreteness and fluctuations in the target (taken to be the larger of the two dipoles) 
wave function||22[|. Although our derivation does not indicate any reason to expect that discreteness 
and/or fluctuations are particularly important, it does indeed show that the application of the 
JIMWLK or KLWMIJ evolution to the target wave function in the dipole-dipole scattering is 
flawed. The reason KLWMIJ evolution fails is obvious. Starting with a dilute single dipole target 
initial stages of the evolution are indeed well described by the KLWMIJ equation. However when 
the density in the target wave function reaches large value j oc 1/g neglecting high density effects 
in the evolution of the wave function is not permissible. This density is parametrically the same 
as that for which the scattering amplitude becomes of order one, and it is therefore also the same 
density at which the effect of the multiple scattering corrections in KLWMIJ evolution becomes 
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significant. This has been recognized in the hterature for a long time, see for example fifth paper 
in 



On the other hand the reason for the failure of JIMWLK is somewhat more subtle. Again 
starting with the dilute target one can initially evolve it with the JIMWLK equation. The multiple 
scattering effects are not important as long as the density is small, and thus the use of JIMWLK 
in the dilute regime is as good as the use of KLWMIJ. When the density is parametrically large 
again the JIMWLK evolution is valid, since the evolution of the wave function is accounted for 
appropriately and the long range multiple scattering corrections dominate at high density. It might 
therefore seem that JIMWLK equation can be used all the way through in this situation. This is 
however not the case. The reason it fails is that there is a range of rapidities in the evolution when 
the density is already not very small but the rate of growth is still large. This happens just before 
the saturation is reached. Since the density in this range of rapidities still grows exponentially, the 
short range multiple scattering effects dominate. Those are not included in JIMWLK evolution, 
and thus the rate of growth of the amplitude is overestimated. Note that if already at the initial 
rapidity the density in the target wave function is large (e.g. for a heavy nucleus) there is no rapidity 
window in which the short range multiple scatterings dominate, and thus JIMWLK evolution is 
valid. 

We close the discussion by stressing that the calculation of the wave function given in the 
present paper is the correct starting point for derivation of the complete evolution equation which 
takes into account all relevant Pomeron loop effects. The validity of such equation will not be 
limited to the process of collision of two small objects, but more interestingly to the situation where 
two colliding objects are large. The use of JIMWLK evolution in this case is not justified since the 
soft gluons produced in the wave function can multiply rescatter on the large target field. 
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